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The superstatistics concept is a useful statistical method to describe inhomogeneous complex sys- 
tems for which a system parameter /3 fluctuates on a large spatio-temporal scale. In this paper we 
analyze a measured time series of wind speed fluctuations and extract the superstatistical distribu- 
tion function directly from the data. We construct suitable Langevin and Fokker-Planck models 
with a position dependent /3-field and show that they reduce to standard type of superstatistics in 
the overdamped limit. 
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The superstatistics concept introduced in 2003 [1] provides a useful tool to describe a large variety of complex 
systems |2h15I|. The basic idea is to characterize the complex system under consideration by a superposition of two 
^ ■ statistics, one corresponding to ordinary statistical mechanics (on a mesoscopic level modeled e.g. by a Langevin 
equation) and the other one corresponding to a slowly varying parameter f3{x, t) of the system which can be, but need 
not to be, an inverse temperature. There are many interesting applications of superstatistical techniques. Recent 
work in this direction includes train delay statistics [l6[ , the distribution of accelerations of test particles in turbulent 
CN . flows and cancer survival statistics [isj . Further applications are described in (Tol - fiH . 

' In this paper, as a working example, we study measured statistics of wind speed fluctuations. For previous work 
(~| . in this direction, see [H, [2^. We will extract the superstatistical distribution function /(/3) out of our data set and 
O ' show that in good approximation it is an inverse gamma distribution. As a theoretical model for these and other data 
, sets, we will study Langevin and Fokker-Planck equations with a varying, position-dependent /3. We will show that in 
^ ■ the overdamped limit these models reduce to standard type of superstatistics. Skewed stationary distributions arise 
I ' , , naturally in these types of models if the overdamped limit is not completely performed. 

d This paper is organized as follows. In section II we recall the basic concepts of superstatistics, illustrated with the 

I example of a linear Langevin equation. In section III we apply our methods to measured wind speed data, extracting 
^ I various relevant distributions and parameters out of the data set. In section IV we extend our theoretical approach 

■ by studying a Kramers-like equation with spatially varying /3. We show that this model reduces to standard type of 
G I superstatistics in the overdamped limit. An example is studied in section V. Then, in section VI we comment on the 

I ' skewness of the generated distributions. The last section gives a short discussion of the results of this paper. 

■ 

^ : 

Q ; II. LINEAR LANGEVIN EQUATION WITH TIME-DEPENDENT P 

I I Let us study as a simple example the following Langevin equation 

> ■ I 

, where 7 and m are constants and g{t) corresponds to Gaussian white noise with unit variance. v{t) is a stochastic 

C\| ' process that can be associated with the velocity of a Brownian particle. If f3{t) = /3 is a constant, then one ends up 

, with the standard Langevin equation which describes the dynamical process of a Brownian particle with mass m in an 

■ environment with constant temperature 1//3. It is well known that the stationary velocity distribution of a Brownian 
' . I ] particle is a Gaussian distribution with zero mean and variance The relaxation time is r = I/7. In superstatistics, 

J> . one considers a Brownian particle moving through an environment with a slowly fluctuating temperature field whose 
changes take place on a typical time scale T such that t « T. When this inequality holds, the local velocity 
distribution of the system can relax to a Gaussian distribution before the next change of f3{t) takes place. As 
such, after a long time, the stationary velocity distribution P{v) of the particle is just a superposition of Gaussian 
distributions weighted with a function 



P{v) ^ J d/3/(/3)y^exp( --m/3«^ 1 . (2) 
'^min 

This /(/?) is the probability density to observe some value of /?. Depending on the properties of /(/?), different 
results for the stationary velocity distribution P{v) will occur e.g., power-laws or stretched exponentials. In 
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[3. [29| a method was introduced to determine /(/?), given an experimental time series. The method extracts the 
main superstatistical parameters out of a given data set and examines the vahdity of the superstatistical model 
assumptions. Depending on the system under study, one can obtain different results for f{/3). We will first briefly 
outline this method to extract from the data. 

The starting point is a discrete time series v containing n data points. Using the definitions proposed in [2|, [29j . 
one can estimate values for the two different time scales r and T. In case that the inequality r/T << 1 holds, one 
proceeds by dividing the time series v in N different time slices of length t with N — \ n/t\ , where [x\ means rounding 
the value of x to the nearest lower integer. Then, one calculates the variance of v in each of these sub-intervals. The 
inverse of the variance is an estimator for the value of /3. As such, one obtains a new series containing N points which 
are denoted as Pi with i = 1, 2, ... TV. Within this assumption, the distribution P{v) is approximated by 




Pi-)-^^}^\H^^-''''' ■ (3) 



When N is large enough, one can replace expression ([3]) by ([2]) with /(/?) being the probability density that the value 
of the inverse variance in a randomly chosen time slice of length T equals /3. Notice that the superstatistical approach 
includes two approximations. In the first step, one assumes the existence of two time scales t and T such that in every 
time slice 'local' equilibrium is reached. Then, P{v) can be approximated by In the second step one assumes 
the existence of a distribution /(/?) replacing the summation in expression Q by an integral. Then, P{v) can be 
approximated by ©. 

III. APPLICATION TO WIND SPEED FLUCTUATIONS 

In the following, as a working example, we apply our method to an experimentally measured time series of the 
horizontal component v{t) of wind speeds recorded af the Lammefjord site at a height of 10m (cup 1) during the 
year 1987 [i^. The incremental distribution of us{t) := v{t + 5) — v{t) exhibits non-Gaussian behavior that can 



be modelled using superstatistics |41| . As an illustration, we calculated the superstatistical approximation of the 
incremental distribution {5 — 8) of the windspeeds at day 191. The measuring frequency is 8Hz. This means that the 
number of data points is rt = 24 x 3600 x 8 ~ 7 x 10^. The two different time scales extracted from the data are t « 4.1 
and T Ki 112, obtained by using similar techniques as in 0, H^. This illustrates that the data set shows clear time 
scale separation, in agreement with the results of [26l |. Knowing the value of T, one can construct the distribution 
/(/?) and try to approximate this histogram with some well-known distributions such as the gamma distribution, the 
lognormal distribution or the inverse gamma distribution. The latter one is given by 

/(/?) = le-^/'^. (4) 

Here a and 9 are parameters. The relevance of the above three distributions was motivated in "2^. The results of our 
calculations are shown in figure [TJ There is an excellent agreement between the histogram extracted from the data and 
the inverse gamma distribution. The right part of the figure also shows the empirical distribution P{u) together with 
the first and second approximation of superstatistics given by eqs. (jS]) and ([2]), respectively. The first approximation 
of superstatistics well models the fat tails of the empirical distribution P{u), see Figure [T] The excellent agreement 
between P{u) and expression with being the inverse gamma distribution, gives further evidence that this 
latter distribution is well able to represent the fluctuations of /?. Summarizing, our simple superstatististical model 
discussed so far is a good first-order approximation of the process of windspeeds increments. 

Figure [5] shows the extracted time scale ratio T/r as a function of 6. Apparently, for small S the time scale separation 
between T and r is less pronounced. The right part of the figure also shows a parameter e as a function of S. This 
e was defined in [29j ) as a measure of quality of the superstatistical approximation. The smaller e, the better the 
superstatistical model assumptions are satisfied for the given time series. 

IV. SUPERSTATISTICS AND OVERDAMPED MOTION 

We will now generalize the Langevin approach to superstatistics. We start from the following set of equations 
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FIG. 1: (Colour online) LEFT: Plot of the empirical distribution /(/3) (dots) extracted from the wind data and the best 
fit to a lognormal distribution (red), gamma distribution (blue) and an inverse gamma distribution (green). RIGHT: Plot 
of the empirical distribution P{u) together with the first superstatistical approximation ^ (black line). Also the second 
approximation ^ is shown, where f(j3) is given by the lognormal distribution (red), the gamma distribution (blue) and the 
inverse gamma distribution (green). 
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FIG. 2: (Colour online) LEFT: Ratio of the time scales T and r as extracted from the data. RIGHT: The parameter e, defined 
in [2^, measures how good a superstatistical model fits a given data set. Ideal superstatistics is corresponding to e = 0. 



dt 



xit) = v{t), 



(5) 



with F{x) = —dV{x)/dx. Here V{x) is a confining potential like, e.g., a harmonic potential, and x is a, position 
variable. Crucial for our approach is that /3 depends on x but not on the velocity v. We are interested in the 
stationary velocity distribution P{v). It is well known that in case I3{x) = /3 is a constant, this distribution becomes 
a Gaussian distribution 



P{v) 




exp 



1 



(6) 



However, we are interested in the general problem where the temperature is position dependent. Then, in order to 
obtain an exact expression for P{v), one has to solve the corresponding stationary Fokker-Planck equation for the 
joint distribution P{v, x), 



= 



d f 1 ,A d 7 1 52 

— \-fv F(x) - V— + —-TT-^-^ 
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P{v,x) 



(7) 
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and integrate out the dependence of the variable x 

/■^max 

P{v) = / dxP{v,x). (8) 

J X 

mm 

The problem with this procedure is that usually it is not possible to obtain an analytical expression for P{v,x). 
However, in order to make the connection with the superstatistical approach outlined in the previous section, it is not 
necessary to have the exact form of P{v, x). It is sufficient to study (0 in the high friction limit, i.e., 7 00. Taking 
this overdamped limit of ([5]) results in a stochastic differential equation in the variable x(t) only: 



= J_F{x) + J—-^g{t) (9) 
at ]l rwy pyx) 

The associated Fokker-Planck equation for the stationary distribution Po{x) becomes 

o = -|:n^)^o(x) + ^^Po(x) (10) 

(we used Ito's interpretation). The lower index o is used to denote the overdamped limit. Equation ([TU)) can be 
solved analytically 

= z-.«.)»p(/..™). (11, 

where Z is & normalization constant. Taking the overdamped limit physically means that the velocity will thermalise 
very quickly. As a consequence, for large values of 7, the solution P(w, x) of the complete Fokker-Planck equation ([7]) 
is approximated by [13] 

P{v, x) « Po{v, x) = Po{x)Po{v\x). (12) 

Here Po{v\x) denotes the conditional distribution of v given the local value of the slow variable x. This is just a 
Gaussian distribution ^ with j3 replaced by I3{x). As a consequence, one obtains the following approximation for 
the distribution of interest: 



=max /-^max I mf3{x) ( 1 



P{v)^J dxPoiv,x)^J dxPoix)^ ^^^ejcpl-^ml3{x)v^j . (13) 

^min ^min 

After a change of variables, one ends up with ([2]) in which the superstatistical distribution f{/3) is associated with 



( dm 

\ dx 



-1 



Poim)- (14) 



This shows that the superstatistical approximation is formally equivalent to taking the high friction limit of a stochastic 
differential equations of type ([S]). The latter model describes the dynamics of a Brownian particle in an inhomogeneous 
heat bath in the presence of a confining potential. 

V. EXAMPLE 

To illustrate the theoretical considerations of the previous section, we study an explicit example here. Assume for 
simplicity that the temperature is linear in the position x. In order to ensure that the temperature remains positive 
we introduce an arbitrary small positive constant a and define 



1 



a 



1 2 



Ixl + a 2 



With this choice one obtains for the velocity distribution in the overdamped limit ([T 



P{v) « J^"^ dxP{x)exp(^~ J dx'^^^m)^ y^^^exp (^-lmf3{xy 
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FIG. 3: (Colour online) LEFT; The difference between the empirical distribution P{u) and the second approximation of 
superstatistics. RIGHT: Kurtosis and skewness of the empirical distribution P{u) as a function of S. 




2tj. exp f -im^(a;)w2 



(16) 

where we absorbed constants in the normalization Z . In the limit a — > 0, this integral simplifies to 

P(.) - ^/f (/fn) ' K„_,/, (V2^h) (17) 

where K^{x) is the modified Bessel function of the second kind. Notice that f{(3) becomes an inverse gamma dis- 
tribution in this limit. This particular example connects the results of sections III and IV because it shows that 
the theoretical model of section IV together with the choice p3|) for the temperature and the potential reproduces 
the observed empirical distributions discussed in section III. However, this is only one possibility. In principle, for 
every choice of /3(x), a possible candidate for the force that results in an inverse gamma distribution for /(/3) can be 
calculated by evaluating the following expression 

Where the accent means differentiating with respect to x. The reason to chose the particular form of T{x) and V{x) 
as in is that this is probably the easiest choice that both fulfills equation ^TE\\ and is physically meaningful. The 
latter is also the reason that we introduced the constant a in the expression for the temperature (|15p . This is to ensure 
that the particle can pass the origin in the overdamped limit. This problem disappears if one takes finite friction into 
account. Therefore, it is appropriate to restrict our analysis to the special case a — >■ 0. 



VI. SKEWNESS 



We showed in section III that the superstatistical approach is well suited to describe the fat tails that are exhibited 
in the wind speed distribution under study. However, it is well known, and clearly visible in figure [TJ that the 
distribution P{u) is slightly skewed (with vanishing first moment). Clearly, the superstatistical distribution ^ is 
symmetric and can not represent the skewness that occurs in the experimental data. Figure [3] (right part) shows the 
skewness and kurtosis as a function of S for the data set of measured velocity increments us{t). It is possible to model 
the kurtosis with the superstatistical approach described up till now, however it is impossible to model the skewness. 

As mentioned before, the complete Fokkcr-Planck equation ([7]) is usually not analytically solvable. However, 
relations between the moments of the distribution P(v,x) can be obtained, using the Fokker-Planck operator L. One 
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has 

(Lz;V)=0 for ^, j ^ 0,1^2, .. .withL := {-,v - -F{x)) — ~ v— ^ — . (19) 

\ TO J OV OX m p[x) o^v 

Applying this with i = 0, j = 1 results in (v) = 0, while depending on the particular choices of F{x) and I3{x) one can 
obtain (v^) ^ 0. This means that the stationary velocity distribution P{v) can have non- vanishing skewness, while at 
the same time it has a vanishing first moment. However notice that the overdamped approximation (fT^ of P{v,x) 
is symmetric in v, regardless of the choice of /3(x). This means that in order to introduce skewness, one has to go 
beyond the overdamped limit. Perturbation expansions of P{v, x) around the overdamped limit can be performed 
|42| . The first order correction is 

P(v,^) = y^^exp (^-iTO/3(x)z;2) {\ + -\ {mfi{xV - 3] ■ (20) 

Notice that the contribution of the first order correction vanishes for with i = 0, 1,2. The third order moment 
becomes 

(.3) ^ r^dxv'Piv,x) = dxp^ix)-^^. (21) 

This shows that the theoretical model introduced in section IV can be used to systematically describe the experimental 
data. The superstatistical approximation (the overdamped limit) is used to model the fat tails (the kurtosis) of the 
empirical distribution. By going beyond the overdamped limit, skewness is naturally introduced in the theoretical 
distribution. 



VII. DISCUSSION 



In previous papers dealing with superstatistics, the superstatistical concept was completely probabilistic in nature. 
One was facing the question what the relevant temperature distributions f{(3) are. In this paper, we have shown 
that the superstatistical model is equivalent to a dynamical process with a position dependent temperature field in 
the overdamped limit. This opens up the possibility to relate the choice of /(/3) to a concrete question with a clear 
physical interpretation, namely what are appropriate choices for the potential and the temperature fields. 

One such example was studied in this paper. If one assumes a temperature field that is linear in the position, in 
combination with basically a harmonic potential, one ends up with an inverse gamma distribution for /(/3). This 
distribution is indeed observed in empirical data of wind speed fluctuations. We emphasize that other choices for 
temperature field and potential are of course possible. It is an interesting topic for further research to examine other 
physically meaningfuU combinations of /3(x) and V{x). One can also try to extract the position dependence of I3{x) 
and V{x) immediately out of empirical data, see e.g. [4^ where such analysis is performed for financial time series. 

Finally, we showed that the theoretical approach to superstatistics presented in this paper naturally leads to skewed 
probability distributions if one goes beyond the overdamped limit. Slightly skewed distributions are observed in several 
scientific fields for various variables, e.g., wind speed fluctuations (see section VI), log returns of prices in the stock 
market [i^ , and velocity increments of turbulent flows [i^ . Apart from the observed skewness in the data, these three 
examples have in common that the symmetric part of the empirical distribution can be well approximated by ^ with 
different expressions for /(/3), see, e.g., the present paper and 29]. In the theoretical model developed in this paper, 
expression (0) is the zeroth-order approximation of the true superstatistical distribution which is obtained by fully 
taking into account the position dependence of I3{x) and V{x). Therefore, the current approach is very promising, 
since one can now go beyond ^ and try to calculate successive higher-order corrections in order to describe all 
the essential features of an empirical distribution at hand. Apart from the aforementioned examples, also dense 
granular flows are interesting experimental systems from a superstatistical point of view, because position dependent 
temperature profiles can be estimated for these kinds of systems. Also, it has been shown that the particle velocity 
distribution shows skewness 
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